Nucleon distribution amplitudes from lattice QCD 
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We calculate low moments of the leading-twist and next-to-leading twist nucleon distribution amplitudes on 
the lattice using two flavors of clover fermions. The results are presented in the MS scheme at a scale of 2 GeV 
and can be immediately applied in phenomenological studies. We find that the deviation of the leading- twist 
nucleon distribution amplitude from its asymptotic form is less pronounced than sometimes claimed in the 
literature. 



Introduction. — Distribution amplitudes [1, 2, 3, 4, 5, 6] 
describe the structure of hadrons in terms of valence quark 
Fock states at small transverse separation and are required in 
the calculation of (semi)exclusive processes. A simple pic- 
ture is obtained at very large values of the momentum transfer. 
For example, the magnetic Sachs form factor of the nucleon 
Gm{Q 2 ) can then be expressed as a convolution of the hard 
scattering kernel h(xi, yi,Q 2 ) and the leading-twist quark dis- 
tribution amplitude in the nucleon ip(xi, Q 2 ) [3]: 



G M (Q 2 ) 
e'2 



' [dx] I ^^(y^Q^hix^y^Q^ix^Q 2 ) 



where [dx] = dxidx 2 da;3i5(l — X^=i x i)> an d — Q 2 i s tne 
squared momentum transfer in the hard scattering process. 
However, in the kinematic region lGeV 2 < Q 2 < 10 GeV 2 , 
which has attracted a lot of interest recently due to the JLAB 
data [7, 8] for Gm, the situation is more complicated. Here 
calculations are possible, e.g., within the light-cone sum rule 



approach [9, 10]. They indicate that higher-twist distribution 
amplitudes become important while higher Fock states do not 
play a significant role. In any case, the distribution amplitudes 
are needed as input. 

Being typical nonperturbative quantities, distribution am- 
plitudes are difficult to compute reliably in a model-indepen- 
dent way. Determinations by QCD sum rules have been at- 
tempted, but suffer from considerable systematic uncertain- 
ties, especially for lower values of Q 2 . As advocated in the 
pioneering work [11], lattice QCD can provide valuable addi- 
tional information. 

In this paper we present an improved and extended lattice 
analysis of the nucleon distribution amplitudes. We find that 
the asymmetry of the leading-twist amplitude is smaller than 
in QCD sum rule calculations, in agreement with phenomeno- 
logical estimates [12, 13], which suggest a less asymmetric 
form. 

General Framework. — In the case of the proton, the 
starting point is the matrix element of a trilocal quark oper- 
ator, 
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where path ordering is implied for the exponentials, a, 6, c are 
the color indices, N the proton spinor and \p) denotes a proton 
state with momentum p. We will consider this matrix element 
for space time separation of the quarks on the light cone Zi = 
a t z {z 2 = 0) and J^i = l - 



In momentum space we have 
V{ Xl )^ Jv{z i -p)\{e^{ix i {z i - P ))^-^- (2) 

i— 1 

with V(xi) = V(xi, X2, x%) and similarly for A(xi) and 
T(xi). The distribution amplitudes V(xi), A(xi) and T(xj) 
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describe the quark distribution inside the proton as functions 
of the longitudinal momentum fractions x- t . The dependence 
on the renormalization scale is suppressed for notational sim- 
plicity. 

We consider the moments of distribution amplitudes, which 
are defined as 



V 



[dx] x\x^x^ V(x 1 ,x 2 ,x 3 ), 



(3) 



with analogous definitions for the other distribution ampli- 
tudes. Using eqs. (1) and (2) one can relate the moments of 
the leading-twist nucleon distribution amplitudes to matrix el- 
ements of the local operators 



x [i n V Ul . . .V^ lb d(0)] c T e abc 



by (0\VP mn (0)\p) = -f N pP l7nn N T V lmn with similar rela- 
tions [11] for the operators _AP lmn and %P lmn corresponding 
to the moments A lmn and T lrnn , respectively. Here the multi- 
indices Ifhn denote the Lorentz structure connected with the 
covariant derivatives on the r.h.s. 

Due to the presence of two Mp-quarks in the proton and 
isospin symmetry, the three different amplitudes can be ex- 
pressed in terms of the single amplitude (f>(xi) with the corre- 
sponding moments 
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(4) 



The normalization constant /jv is defined by the choice 
^ooo _ i -pjjg momen t s f the combination ip(xi) = 

V(xi) — A(xi), usually used in sum rule calculations, can 
easily be obtained as ip lrnn — 2<fi lmn — ^ nml _ i n the numer- 
ical calculation, however, we prefer the combination (fr lmn as 
the corresponding statistical errors are smaller by a factor of 
about 3. Note that momentum conservation implies 



In particular we have 



1 = 
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In the limit of Q 2 — > oo one gets <f>(xi) = 120xiX2a;3 [4] and 
the moments <j) lmn are known exactly: 100 = <f w = 001 = 
1/3, </> 200 = (j> 020 = </> 002 = 1/7 and 011 = </> 101 = </> 110 = 
2/21. Thus asymmetries of the type 100 — 010 are important 
quantities at low energies as they describe the deviation from 
the asymptotic case. 

In the case of the next-to-leading twist distribution ampli- 
tudes we restrict ourselves to operators without derivatives, 





/3 = 5.29 


P = 5.40 


fw ■ 10 3 [GcV 2 l 


2.984(60)(157)(65) 


3.144(61)(29)(54) 


-Ai ■ 10 3 [GeV 2 ] 


39.69(76)(259)(124) 


38.72(70) (43) (106) 


A 2 ■ 10 3 [GeV 2 ] 


78.70(155)(562)(245) 


76.23(139)(84)(207) 


0100 


0.3549(11)(61)(2) 


0.3638(11)(68)(3) 


^010 


0.3100(10)(73)(1) 


0.3023(10)(42)(5) 


0001 


0.3351(9)(11)(2) 


0.3339(9)(26)(2) 


0100 _ 0001 


0.0199(23) (46) (4) 


0.0300(23)(93)(1) 


0001 _ 0010 


0.0251(16)(84)(3) 


0.0313(17)(12)(7) 


<A 011 


0.0863(23) (97) (74) 


0.0724(18)(82)(70) 


0101 


0.1135(23) (3) (33) 


0.1136(17)(32)(21) 


<A 110 


0.0953(21)(58)(31) 


0.0937(16)(3)(38) 


0200 


0.1508(38)(213)(64) 


0.1629(28)(7)(68) 


0020 


0.1207(32)(43)(56) 


0.1289(27)(37)(51) 


0002 


0.1385(36)(47)(64) 


0.1488(32)(77)(73) 


0110 _ 0011 


0.0075(33) (69) (44) 


0.0211(27)(78)(32) 


0101 _ 0110 


0.0172(29)(82)(57) 


0.0204(21)(134)(50) 


0200 _ 0020 


0.0335(43)(26)(78) 


0.0321(33)(69)(55) 


0002 _ 0020 


0.0170(36)(8)(56) 


0.0193(24)(32)(42) 



TABLE I: Moments and asymmetries in the MS scheme at 2 GeV 
for 4> lmn = {V lmn - A lmn + 2T lnm )/3. The first error is statis- 
tical, the second (third) error represents the uncertainty due to the 
chiral extrapolation (renormalization). The systematic errors should 
be considered with due caution, see the text for their determination. 



i.e., to the lowest moments. Thus the problem is simplified 
greatly since the Lorentz decomposition of the relevant ma- 
trix element involves only two additional constants Ai and A2 
[14]. They describe the coupling to the proton of two indepen- 
dent proton interpolating fields used in QCD sum rules. One 
of the operators, C T , was suggested in [15] and the other, Ai T , 
in [16]: 

C T (0) = e abc [u aT (0)CYu b (0)\ x ( 75 7prf c (0))r, 
MM = e abc [u aT (0)Ca^u»(0)] x ( l5 a^d c (0)) T . 

Their matrix elements are given by 



(0|£ T (0)|j>) = X imN N T , 
(0|M T (0)|p) = X 2 m N N T . 



(7) 

(8) 



Due to Fierz identities we have 

8 
m N 



(2Ai + X 2 )N = —(0\e abc (u aT Cd b ) l5 u c \p), (9) 



which vanishes in the nonrelativistic limit. 

Computation. — The required matrix elements between 
the vacuum and the proton state are extracted from two-point 
correlation functions with the investigated local operators at 
the sink and a smeared interpolating operator for the proton at 
the source. In addition one needs the usual proton correlator 
with both source and sink smeared. We have evaluated these 
two-point functions on gauge field configurations generated 
by the QCDSF/DIK collaborations with the standard Wilson 
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gauge action and two flavors of nonperturbatively improved 
Wilson fermions (clover fermions). The gauge couplings used 
are j3 = 5.29 and (3 = 5.40 corresponding to lattice spacings 
a ps 0.075 fm and a ps 0.067 fm via a Sommer parameter 
of ro = 0.467 fm [17, 18]. Our smallest pion masses are 
380 MeV {(3 = 5.29) and 420 MeV (f3 = 5.40), while the 
spatial lattice sizes L are such that m^L > 3.7. 

Due to the discretization of space-time, the mixing pattern 
of the operators on the lattice is more complicated than in the 
continuum. It is determined by the transformation behavior 
of the operators under the (spinorial) symmetry group of our 
hypercubic lattice. As operators belonging to inequivalent ir- 
reducible representations cannot mix, we derive our operators 
from the irreducibly transforming multiplets of three-quark 
operators constructed in [ ; ] in order to reduce the amount 
of mixing to a minimum. These irreducible multiplets con- 
stitute also the basis for the renormalization of our operators, 
which is performed nonperturbatively. To this end we contract 
our three-quark operators with three quark sources, amputate 
the external legs from the resulting four-point functions and 
impose an RI'-MOM-like renormalization condition. Finally 
we use continuum perturbation theory and the renormaliza- 
tion group to convert the results to the MS scheme at a scale of 
2 GcV. We estimate the corresponding uncertainty by varying 
the scale at which our renormalization condition is imposed 
between 10 GeV 2 and 40 GeV 2 . In this procedure, the mixing 
with "total derivatives" is automatically taken into account. 

In the case of the moments considered in this work we can 
avoid the particularly nasty mixing with lower-dimensional 
operators completely. Note that the operators y^' mn , j[pimn 
and T T plmn with different multi-indices plfhn but the same 
Iran are related to the same moments V lmn , A lmn and T lmn , 
and we make use of this fact not only in order to minimize 
the mixing problems but also in order to reduce the statistical 
noise by considering suitable linear combinations. 

For the operators without derivatives, i.e., the matrix ele- 
ments Ai, A2 and Jm, we have performed a joint fit of all 
contributing correlators to obtain the values at the simulated 
quark masses. As these are larger than the physical masses a 
chiral extrapolation to the physical point is required in the end. 
To the best of our knowledge there are no results from chi- 
ral perturbation theory to guide this extrapolation. Therefore 
we have adopted a more phenomenological approach aiming 
at linear (in to 2 ) fits to our data. It turns out that the ratios 
In / m % an d ^i/ m N are particularly well suited for this pur- 
pose (see Fig. 1 (upper panel) for an example). Moreover, 
/jv /jtijv is dimensionless and hence not affected by any un- 
certainty in the determination of the lattice spacing. In order 
to estimate the systematic error due to our linear extrapola- 
tion, we also consider a chiral extrapolation including a term 
quadratic in to 2 and take the difference as the systematic error. 
The results in the MS scheme at a scale of 2 GeV are given in 
Table I. Note that 2Ai ps — A2, a relation that is expected to 
hold in the nonrelativistic limit (cf., eq. (9)). 

For the higher moments one can proceed in the same way 
and the constraint (6) is satisfied very well. However, the 
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FIG. 1 : Linear chiral extrapolation of bare lattice results for /jv / m 2 N 
(upper panel) and the asymmetry R 200 — R 020 (lower panel) with 
one-cr error band. 





LAT 


QCDSR 


BLW 


BK 




0.304 


0.248 


0.303 


0.311 


2A oio 


0.091 


0.303 


0.116 


0.064 



TABLE II: Comparison of our results (LAT) to selected sum rule re- 
sults [20] (QCDSR) and the phenomenological estimates [ 1 2] (BLW) 
and [13] (BK) at the scale 2 GeV. 



statistical errors in this approach are too large to allow an 
accurate determination of the (particularly interesting) asym- 
metries. We achieved smaller errors by calculating the ratios 
R imn = where Si = loo + olo + </> 001 forZ + m + 

n = 1, and S 2 = 2(0 011 + 101 + cf> 110 ) + </> 200 + </> 020 + 002 
for I + ra + n = 2. These ratios can be extrapolated linearly 
to the physical masses. An example is shown in Fig. 1 (lower 
panel) for the case of the asymmetry R 200 — i? 020 . Requiring 
that the constraint (6) be satisfied for the renormalized values 
we can finally extract the moments from the ratios. 

Discussion and Conclusions. — Since we only have re- 
sults at two different lattice spacings, we are unable to extrap- 
olate our results to the continuum limit. However, we find 
that the results obtained at (3 = 5.29 and (3 = 5.40 are com- 
patible within errors. Hence we take the data from our finer 
lattice ((3 = 5.40) as our final numbers. The values for (j) lmn 
imply that (p wo = 0.394, ip ow = 0.302 and ip mi = 0.304. 
These moments can be interpreted as the fraction of momen- 
tum carried by the corresponding quarks [5, 6]. Thus we find 
that the largest fraction of the proton longitudinal momentum 
is carried by one u/9-quark with spin aligned with the proton 
spin. However, this asymmetry is not as strong as found in the 
QCD sum rule calculation. Our results for the first moments 
are close to phenomenological estimates [12, 13], cf. Table II. 
On the other hand, our results for ip 011 , (p 101 and ip 110 are 
similar to the sum rule values, while the asymmetries in the 
moments tp 200 , ip 020 and (p 002 are clearly smaller. 

Let us now expand the distribution amplitude in terms of 
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FIG. 3: The model distribution amplitude (p(xi, X2, xz) for X3 
0.5 as a function of x\ with statistical errors. 



FIG. 2: Barycentric contour plot of the leading-twist distribution 
amplitude <p(xi, X2, £3) at fi = /10 = 2 GeV as obtained from the 
moments presented in Table I. The lines of constant xi, X2 and £3 
are parallel to the sides of the triangle labelled by x%, 2:3 and x\, 
respectively. 



polynomials P n to order N 
matrix is diagonal [21, 22]: 



2 chosen such that the mixing 



N 



ip{Xi,(J,) = l20xiX2X 3 y~] C n (fJ, )Pn(Xi) ( 



n=0 



Calculating the coefficients c n (no) from an independent sub- 
set of the moments <fi lmn (fio = 2 GcV), we obtain a 
model function for the distribution amplitude presented in 
Fig. 2. While the (totally symmetric) asymptotic amplitude 
120 X1X2X3 has a maximum for x\ = x-2 = X3 = 1/3, inclu- 
sion of the first moments (i.e., choosing N = 1) moves this 
maximum to x\ ~ 0.46, X2 ~ 0.27, £3 « 0.27 giving the 
first quark substantially more momentum than the others. The 
second moments then turn this single maximum into the two 
local maxima in Fig. 2. The approximate symmetry in X2 and 
X3 is due to the approximate symmetry (p lmn « ^inm Q f our 
results. It is also seen in QCD sum rule calculations as well 
as in several models such as BLW and BK and may indicate 
the formation of a diquark system. To illustrate the statistical 
uncertainty we show in Fig. 3 the profile of ip at 23 = 0.5 and 
the corresponding error band. Note that higher-order polyno- 
mials have been disregarded in this model and thus Figs. 2 and 
3 should be interpreted with due caution. 

In order to establish a link to observable quantities we are 
calculating nucleon form factors using light cone sum rules. 
As our moments are close to the phenomenological estimates 
shown in Table II we can expect reasonable agreement with 
the experimental data. It is, however, to be stressed that these 
calculations still involve some model dependence, while the 
moments presented in this work were obtained from first prin- 
ciples. 
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